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ABSTRACT 

We study a simple aggregation model that mimics the clustering of traffic on a one-lane 
roadway. In this model, each "car" moves ballistically at its initial velocity until it overtakes 
the preceding car or cluster. After this encounter, the incident car assumes the velocity of 
the cluster which it has just joined. The properties of the initial distribution of velocities in 
the small velocity limit control the long-time properties of the aggregation process. For an 
initial velocity distribution with a power-law tail at small velocities, Pq(v) ~ v M as v — > 0, 
a simple scaling argument shows that the average cluster size grows as n ~ 

£ ( M +l)/(M+2) 

and that the average velocity decays as v ~ t -1 /(^+ 2 ) as t — > oo. We derive an analytical 
solution for the survival probability of a single car and an asymptotically exact expression 
for the joint mass- velocity distribution function. We also consider the properties of spatially 
heterogeneous traffic and the kinetics of traffic clustering in the presence of an input of 
cars. 

P. A. C. S. Numbers: 68.70+w, 03.20.+i, 05.20.Dd, 05.40.+j. 
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I. Introduction 



A variety of approaches have been applied to describe the collective properties of 
traffic flows [1]. For example, to mimic congested traffic flow in two dimensions, cellular 
automation models have been proposed [2,3]. Asymmetric hopping processes have also 
been applied to model traffic flow on a one-dimensional road [4-6]. When the number of 
cars is large, traffic flows can be modelled phenomenologically in terms of a one-dimensional 
compressible gas [7-9]. Such an approach predicts the appearance of shock waves, where 
hydrodynamic quantities, such as the average density and velocity, become discontinuous. 
However, the hydrodynamic approach does not naturally describe the behavior of traffic 
flows in the low-density limit where there are large heterogeneities in traffic density. For 
this situation, a microscopic model may provide a more appropriate description. 

In this article, we introduce a ballistic aggregation process to model the kinetics of 
clustering in one-dimensional traffic flows. Our approach is inspired, in part, by the recent 
interesting results that have been obtained for a variety of reaction processes which involve 
ballistic particles including: ballistic agglomeration, Ai+Aj — > A^ + j, with momentum con- 
serving collisions [10,11]; ballistic annihilation, A + A — > 0, [12,13]; and several nucleation 
and ballistic growth processes [14-16]. In our model, cars move ballistically in a one direc- 
tion, say to the right, according to an initial velocity distribution. Clusters form whenever 
a faster car overtakes a slower car or cluster. The overtaking car then assumes the velocity 
of the lead car in the cluster. This model is an idealized description for one-lane traffic 
flow. While there are obvious shortcomings in our model, it is exactly soluble and permits 
a thorough understanding of the kinetics of the aggregation process. 

This paper is organized as follows. In section II, we present the model and postulate 
the scaling behavior for the velocity and the concentration of the clusters. This approach 
makes use of the statistical properties of the minimal random variable within a large sam- 
ple. In section III, we investigate the distribution of cluster velocities. For this distribution, 
the cluster size is irrelevant and this feature allows us to consider a simpler "coalescence 
only" model. For this reduced problem, the velocity distribution is obtained exactly in 
terms of the initial distribution of car velocities and then evaluated for general continuous 
distributions. Building on these results, the general clustering process is solved in section 
IV and an asymptotically exact expression for the joint cluster mass-velocity distribution 
is obtained. In section V, we present a formal solution for the velocity distribution function 
for an inhomogeneous initial distribution of particles. We examine the temporal behavior 
that arises for a simple step function initial spatial distribution. In section VI, we investi- 
gate another generalization of the model to the situation with a spatially and temporally 



2 



homogeneous input of cars. Depending on the functional form of the input velocity dis- 
tribution in the low-velocity limit, the input can give rise to a steady state or to a system 
which continues to evolve indefinitely. We give our conclusions in section VII. The details 
of specific calculations are given in the Appendices. 

II. Scaling Analysis 

We consider an idealized one-dimensional traffic flow in which the size of each car is 
zero. This is appropriate for describing clustered traffic in the low-density limit, a situation 
which is often encountered on rural secondary roads. In the following, we will refer to such 
sizeless cars as particles. We consider the initial condition when there are only isolated 
particles ("monomers") in the system with a random spatial distribution of density c . 
The initial velocity distribution, P(v, t = 0), can generally be written in the scaling form 

P(v,t = 0) = ^P o (-) v>0, (1) 
vo \v J 

with / °° Pq (z)dz = 1. Here we have tacitly subtracted the assumed finite value of the 
velocity of the slowest car from all velocities. In what follows, it is often convenient to 
introduce the dimensionless density c/cq — > c, velocity v/vq — » v and time c^v^t — > t. This 
yields a rescaled initial concentration which is equal to unity. 

In our model, particles move at their initial velocities and whenever a particle overtakes 
a cluster aggregation occurs. The aggregation rule is simply that two colliding clusters form 
a new cluster with a velocity equal to the smaller of the two incident cluster velocities and 
with a mass equal to the sum of the two cluster masses (Figure 1). If we denote a cluster 
of mass m and velocity v by A m „, the process is described by the reaction scheme 

^mi,»i "l~ ^tri2, V2 > J 4mi+m 2 ,mm{wi,t)2}- (2) 

We present now a simple argument, based on the statistics of extremes, to predict the 
asymptotic time dependence of the typical cluster mass m and typical cluster velocity v 
at time t. Since the typical distance, /, between clusters grows with time as / ~ vt, the 
typical number of particles in a cluster is proportional to this distance, yielding m ~ / ~ vt. 
To find the typical velocity, one has to relate the mass of a cluster to its velocity. Such 
a relation may be found exactly for an auxiliary "one-sided" problem in which particles 
are placed with a fixed density c to the left of a given "target" particle which moves at 
velocity v, and no particles are placed to the right. Eventually, this target particle will form 
a cluster that includes all consecutive particles to its left whose initial velocities are larger 



3 



than v. The probability that there are exactly k such particles is equal to -P--P+, with 
P+(v) (P_(t>)) defined as the probability that a particle has a velocity larger (smaller) 
than v, i. e., P+(v) = P (v')dv'. Therefore, the average number of particles in the 
cluster that ultimately forms is given by 

oo 

(m(v)) = J2kP-Pl = P + /P- (3) 
1 

Let us now assume a power-law behavior of the initial velocity distribution for small 
velocities, 

P (v) ~ ai> M , v < 1 (4) 
with fj, > — 1 for normalizability. Imposing this power-law form in Eq. (3) yields 

for sufficiently low velocities. For a particle which moves with the typical velocity, it is 
reasonable to expect that this result for the "one-sided" problem gives the correct behavior 
for the full "two-sided" problem. If we combine Eq. (5) with our previous estimate m ~ vt, 
we find the following asymptotic relations, 

_ n + l 

m ~ t a with a ~ .. _i_ o 

• (6) 

v ~ t~P with a _ _ 

Since the mass is conserved in the aggregation process, the typical cluster mass and 
the concentration of clusters c are related by c ~ 1/m ~ t~ a . Notice that in the limit 
H — > oo, the mass grows linearly with time. In contrast, when /z — > — 1, the mass is roughly 
constant, since the velocity distribution becomes effectively unimodal and collisions are 
exceedingly rare. This qualitative dependence on the form of the initial velocity distri- 
bution is reminiscent of the ballistic annihilation process [13], where ballistically moving 
particles annihilate upon collision. In both processes, one finds that the fundamental ex- 
ponents are related by a + f3 = 1 as a consequence of the relation c ~ 1/vt. Moreover, for 
the two processes the decay exponents have similar functional dependences on the form of 
the initial velocity distribution. However, the values of the decay exponents are different: 
for example, for a uniform distribution (corresponding to fi = 0) one finds a = 1/2 for 
the traffic model, while a = 0.76 is obtained in simulations of the annihilation process. 
Additionally, despite the qualitative similarity between these two models for continuous 
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velocity distributions, different behaviors occur when the velocities are discrete. For such 
discrete distributions, the concentration typically decays algebraically in time for ballistic 
annihilation [13], while the concentration decays exponentially in time for the traffic model. 

Since both the typical mass and the velocity scale as power laws in time, the probability 
of finding a cluster of mass m and velocity v, P m (v, t), is expected to evolve toward a scaling 
distribution. Taking into account mass conservation, dv'^2 rn mP rn (v,t) = const., we 
postulate the scaling form 

p m M)~^- 2 «$(M,n (7) 

for the joint distribution, where the scaled mass, M, and scaled velocity, V, are defined by 

M = m/t a and V = vt 13 . (8) 

Note that while the mass m is a discrete variable, the rescaled mass M is continuous. 

Once the joint mass- velocity distribution is found, the single- variable mass and velocity 
distributions can be found by suitable integrations over the subsidiary variable. Thus the 
velocity distribution, P(v,t) = ^2 m P m (v,t), should have the scaling form 

P{v,t)^t p - a %l){V), (9a) 

with ip(V) = dM<$>(M,V), while the cluster- mass distribution, P m (t) = f Q P m (v,t), 
should have the scaling form 

P m (t)~t- 2a 4>(M), (96) 

with 0(M) = J °° dV$(M, V). 

III. The Car Survival Probability 

As a preliminary step in obtaining a full solution for traffic clustering, consider the 
velocity distribution P(y,t). For this quantity, we can ignore the masses of each cluster 
and focus only on the survival probability of a given car. Thus the evolution of the velocity 
distribution is governed by the "derived" coalescence process 

+ Ay 2 > A Ia [ Q ^ Vl y 2 y. (10) 

In the coalescence process, the density of particles with velocity v is identical to P(y,t), 
the velocity distribution of clusters in the full traffic aggregation model defined by Eq. (2). 

Let S(v,t) be the survival probability of particles of velocity v at time t. Here "sur- 
vival" means a car does not overtake any traffic, but an overtaken car is defined to have 
survived. Then the velocity distribution function is given by 

P(v,t) = P (v)S(v,t). (11) 
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The survival probability S(v, t) can be found by considering the possible collisions of a 
particle with initial position and velocity (x, v) with slower particles whose initial positions 
are to the right of x. A collision between the initial particle with co-ordinates (x, v) and 
a slower w'-particle does not occur up to time t if the interval [x, x + (v — v')t] does not 
include the slower particle. For an initial velocity distribution, Po(v), and a Poissonian 
initial spatial distribution, the probability that there is no particle with velocity between 
v' and v' + dv' in the interval [x, x + (v — v')t] is 

exp[-dv'P (v')(v-v')t]. (12) 

For a particle to survive to time t, this exclusion probability should be taken into 
account for every v' < v. To verify this, let us assume otherwise and derive a contradiction. 
Thus consider a particle with initial data (x, v) that has maintained its original velocity 
to time t. In addition, assume that a slower w'-particle is initially present in the above 
exclusion zone, i. e., Ax(0) < Av(0)t. Here Ax(t) is the distance between the two particles 
and Av(t) the relative velocity at time t. Since the velocity v' can only decrease over time 
due to collisions, one has Av(t) > Av(0). Consequently, at time t, the separation between 
the two particles, Ax(t) = Ax(0) - f* Av(t')dt' < Ax(0)-Av(0)t < 0. Thus the i>-particle 
does not survive, in contradiction with the original assumption. 

Hence, the survival probability is simply a product of the exponential factors of 
Eq. (12) for all v' ', with v 1 < v. Evaluating this product gives the survival probability 



S(v, t) = exp 



-t [ dv'(v -v')P Q [v') 
Jo 



(13) 



and combining with Eq. (13) yields the velocity distribution, 



P(v,t) = P (v)exp 



-t [ dv'(v-v')P (v') 
Jo 



(14) 



This is valid for an arbitrary initial velocity distribution P (v ); the only source of stochas- 
ticity arises from the initial conditions. For discrete initial distributions it is seen from 
Eq. (14) that the approach to the final concentration is exponential in time. Thus, we 
focus only on the more interesting continuous initial velocity distributions. 

For the power-law initial velocity distribution Po(^) — av ^ for t> <g; 1, a direct calcula- 
tion shows that the long-time velocity distribution approaches a form that is independent 
of the details of the large- velocity tail of the initial distribution, 

P(v, t) ~ av^ exp [-btv^ 2 ] , (15) 
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with b = a/(fji + + 2). This expression can be written in the scaling form (9a) with 
the scaling function 

^(V) = aV» exp [-bV^+ 2 ] . (16) 

From this solution we see that the velocity distribution maintains the original power-law 
form for small velocities. The exact solution also validates the scaling assumption that the 
asymptotic decay as well as the shape of the limiting distribution are determined solely by 
the low-velocity tail of the initial distribution which, in turn, is governed by the exponent 
fi. 

From Eq. (15), it is straightforward to compute the total concentration and the average 
cluster velocity, 

r 00 if 00 

C(t)= dvP(v,t), (v(t)) = —r dwP(v,t). (17) 

Jo c(t) Jo 

This gives 

c(t) ~ (iJ + l)bT(a)(bt)- a , (18) 

and 

respectively. These expressions confirm the scaling laws suggested in Eq. (6). 

Interestingly, the exact solution of Eq. (14) satisfies the following Boltzmann-like 
integro-differential equation, 

9P(v,t)_ , ^ „.„ D ,„, n , (2Q) 



dt 



= -P(v,t) [ dv'(v -v')P(v',0). 
Jo 



This equation suggests that the loss of v-particles due to collisions with slower i/-particles 
occurs at a rate proportional to the relative velocity, (v—v 1 ). Moreover, the pair correlation 
function factorizes into a product of single-particle velocity distributions, P(v, v', t) = 
P(v, t)P(v', 0) but with different time arguments for the two factors. In contrast, in the 
conventional Boltzmann equation, the decomposition would involve the same argument for 
each velocity distribution. Thus, the exact Eq. (20) quantitatively indicates the degree of 
approximation of the mean-field Boltzmann equation. 

IV. The Full Probability Distribution 

We now solve for the joint mass-velocity distribution function for the general traffic 
model. To obtain P m (v, t), the density of clusters of mass m and velocity v, it is useful to 
introduce the cumulative distribution, Q m (v,t), the distribution of clusters of velocity v 



7 



and mass greater than or equal to m. Once the latter distribution is known, P m (v,t) can 
be obtained by 

Pm(v,t) = Q m (v,t)-Q m+1 (v,t). (21) 

Notice that the density of clusters of mass greater or equal to one is equal to the total 
cluster density, Qi(v,t) = P(v,t). 

Consider a cluster of velocity v which contains at least m particles. Let us number the 
consecutive particles in a cluster from right to left by the index i and denote the rightmost 
particle as i = 0. Denote the initial distance between the i th and (i — l) th particle as x^, 
as illustrated in Fig. 2. We first solve for Q m =2(v,t) and then generalize to any m. Since 
Q 2 (v, t) is the probability that a cluster of velocity v has at least two particles at time t, it 
is equal to the product of the probability that the particle i = has survived up to time t, 
P(v, t), and the probability that the cluster i = 1 (whose mass may be larger than unity) 
collides with the particle i = prior to time t. 

For this collision to occur, the collision partner from the left (i = 1) must have a 
velocity larger than v and the interval x\ < (v\ — v)t must be free of other clusters. The 
probability for this composite event is simply the product of each individual event. Since 
an interval of length x\ is empty with probability exp(— x\), the collision probability is 



Q 2 (t) = P(v,t) J dviPo(vi) J diciexp(-xi). (22) 

xi<(vi— v)t 

The fact that the ^i-particle cannot be slowed down by any other particle before colliding 
with the i>-particle is crucial in obtaining the solution. 

To derive Q m (v, t) for general m, the joint velocity distance distribution Po(vi) exp(— Xi) 
is integrated over the position and velocity of the i th particle for i = 1, . . . ,m— 1. To ensure 
a collision, all m — 1 particles have to move faster than the lead particle and the distance 
of the i th particle from the lead particle must obey x\ + • • • + < (vi — v)t. Imposing these 
constraints on the integration over the velocity and distance of the n — 1 trailing particles 
yields the formal exact expression for the cumulative mass-velocity distribution, 

(.),„{ c.l) Pie. I) ] i / dviP (vi) / dxiexp(-Xi). (23) 



i=i 



xiH \-Xi<(vi—v)t 



For the initial velocity distribution given by Eq. (4), we find the following asymptotic 
behavior (see Appendix A) 

Q m (v, t) ~ tP-W exp [-b(V + MY +2 ] , (24) 
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in terms of the scaling variables M = m/t a and V = vt 13 . 

In the long time limit, the joint mass-velocity distribution, P m (v,t) = Q m (v,t) — 
Q m +i(v,t) can be approximated by P m ~ —dQ m /dm. Performing the differentiation 
gives, 

P m (v, t) ~ t p ~ 2a a b(/i + 2)V»{V + MY +1 exp [-b(V + M)^+ 2 ] , (25) 

which has been explicitly written in the asymptotic scaling form of Eq. (7). This result 
provides a complete description of the traffic aggregation process. It may be considered as 
the ballistic counterpart of the well-known result [17] for diffusion-controlled aggregation 
in one dimension. 

For arbitrary \i we are unable to evaluate the integral over the velocity and obtain the 
explicit mass distribution. However, for the particular case of a uniform initial velocity 
distribution, fj, = 0, it is straightforward to show that 



</>(M) = aexp(-aM 2 /2) ^ = 0. 



(26) 



Another amusing feature of the joint distribution function <fr(M, V) for \x = is the sym- 
metry with respect to the variables V and M. Thus the cluster mass distribution, 0(M), 
and the cluster velocity distribution, i(j(V), are identical Gaussian functions. 

Generally, we are able to extract only asymptotic behavior from Eq. (23). However, 
in the special case of exponential distribution, Pq(v) = e~ v , one can perform all the 
integrations and obtain an explicit solution as detailed in Appendix B. 

V. Clustering in Heterogeneous Traffic Flow 

The above approach can be generalized to the case of a spatially heterogeneous initial 
velocity distribution, Pq(x,v). For simplicity, we ignore the masses of clusters and limit 
ourselves to studying the velocity distribution. This time and space dependent velocity 
distribution, P(x,v,t), may be found by a straightforward generalization of the approach 
developed in section III for the spatially homogeneous case. The resulting expression for 
P(x, v, t) reads 



pV rX—v t 

P(x, v , t) = Pq{x — vt, v) exp — / dv' I dx'Po(x\v') 

JO Jx-vt 



(27) 



As an illustrative example of the effects of a heterogeneous initial particle distribution, 
consider the one-sided distribution in which particles are placed with a fixed density to 
the left of the origin and there are no particles to the right. Thus Pq(x, v) = 9(—x)Pq{v), 
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with 6 the Heaviside step function. For this initial distribution, Eq. (27) yields 

-x/t 



P(x, v, t) = 6(vt — x)Pq(v) exp 



pv p x / z 

-t dv'(v - v')Po(v') - (vt - x) / dv'P (v') 

J x/t JO 

"(28) 

In the long-time limit, the average velocity decays as and hence the front propa- 
gates as x = vt ~ t a . Since the velocity and the position of the front scale as power laws in 
time, Eq. (28) can be expected to have a scaling form. Indeed, by introducing the scaled 
variables X = x/t 01 and V = vt 13 , one can recast Eq. (28) into the scaling form 

P(x,v,t) ~^- a *(X,V), (29) 

with the scaling function 

*(X, V) = 6{V - X)aV exp [-b{V+ 2 - X» +2 )] . (30) 

When X = 0, this scaling function coincides with Eq. (16), the velocity distribution 
for the homogeneous case, *(0, V) = ip(V). Notice also that the density of clusters at 
scaled position X, c(X,t), equals 

pOO I poo 

c(X,t) = c(t)J dV^{X,V) J dViP(V), (31) 

with c(t) given by Eq. (18). In the large X limit, Eqs. (30) and (31) yield 

c(X,t) ~ for X>1. (32) 

X 

Consider now the total number of clusters that infiltrate the initially empty positive 
half-line, N(t) = t a J °° c(X, t) dX. The asymptotic behavior of this quantity is actually 
determined by the finite upper cutoff of the integral, which in turn is given by the position 
of the rightmost particle. For such particles the velocity is of order unity and hence 
Supper = £maxA a ~ vt/t a ~ t& . Therefore 

N(t)~ ^~^ dX ~ al °s(t)- (33) 

Thus the number of clusters entering the empty half-line grows only logarithmically with 
time for arbitrary initial velocity distributions. The only dependence on Pq(v) in Eq. (33) 
is the prefactor a = {ji + l)/(m + 2). 
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VI. Clustering in Traffic Flow with Input 

In this section, we investigate traffic clustering when there is a spatially uniform 
input of cars. This generalization is motivated by real traffic where cars may enter and 
exit a roadway. For the specific case of a spatially homogeneous input of cars we can 
determine the velocity distribution using techniques similar to those employed for the 
traffic coalescence model with no input. 

Denote by R(v, t) the input rate of particles with velocity v at time t per unit length. 
The velocity distribution function for this system, P(v, £), can be expressed as a convolution 
of the flux and the the probability that a particle which was injected at time t' maintains 
its velocity v up to time t in the presence of the input, Sj(v, t, t'), 

P(v,t)= [ dt' #(v,£')SjM,0- ( 34 ) 
Jo 

In writing this expression, we have assumed that the system is initially empty. A particle 
which was injected at time t' will survive until time t if it avoids collisions with all slower 
particles which were present in the system at time as well as avoids collisions with all 
particles which are injected at later times t" > t'. The probability for this composite event 
is simply the product of the probabilities of each event, 



Si(v, t, t') = exp 



-jf dv'{v-v')^P(v',t'){t-t') + £ dt"R(v',t"){t-t")^ 



(35) 



Here the first factor, obtained from Eq. (13) by replacing Pq(v') with P(y',t'), yields 
the probability of avoiding collisions with particles injected prior to time t' . The second 
factor represents the product of exclusions of the type off Eq. (12) for times larger than 
t' and velocities smaller than v' . This factor accounts for the probability that there are 
no collisions with particles injected at times t",t" > t' . Note that the kernel of the 
second factor involves the input rate R at time t", which plays the role of the initial 
velocity distribution at this instant of time. Substituting Eq. (34) into Eq. (35) gives a 
nonlinear integral equation that describes the kinetics of traffic clustering in the presence of 
homogeneous particle input. In the following, we will assume that the flux is constant both 
in time and space, R(v,t) = Pq(v) with f Q dvPo(v) = 1, so that the governing integral 
equation is 



pt r -i rv pv 

P(v,t) = P (v) dt'exp -~t' 2 dv'(v-v')P (v)-t' dv'{v - v')P(v', t - t') 
Jo L 2 J Jq 



(36) 
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In parallel with the case of no input, consider again an initial velocity distribution 
with a power-law small- velocity tail, Po(v) ~ u M . Furthermore, let us assume that the 
concentration and velocity continue to vary as c ~ t~ a and v ~ respectively, and that 
the velocity distribution continues to have the scaling form of Eq. (9), P(v, t) ~ t^~ a ip(V). 
Substituting these into Eq. (36), one can extract consistency conditions for the exponents 
a and (3. For example, the powers of time on both sides of Eq. (36) should be equal - hence 
(3 — a = — Pfj, + 1. Furthermore, both terms in the exponential in the right-hand side of 
Eq. (36) cannot depend on time explicitly - hence (3{jJL + 2) = 2 and a + /3 = 1, respectively. 
These conditions yield a = + 2) and f3 = 2/(/j, + 2). Notice, however, that when [i is 
positive the exponent a is also positive and therefore the concentration, c ~ t _a , decays to 
zero. This is in obvious contradiction with the nature of the problem: with a constant flux 
the concentration may grow indefinitely or a steady state concentration may be reached. 
Thus, one can expect that the above description holds only for \i < 0. For positive \i, 
we anticipate that the system reaches a steady state with a constant concentration and a 
typical cluster mass n ~ t. At the transition fx = 0, a logarithmic temporal dependence is 
anticipated to occur. 

While we cannot confirm the above picture rigorously, we can provide heuristic justi- 
fication. First assume that the velocity distribution evolves towards a steady state Poo(t'). 
Then as t — > oo Eq. (36) becomes 



poo r -i pv pv 

P oo (v)=P (v) cft'exp --t' 2 dv'(v-v')P (v)-t' dv'{v - u')Poo(v') 
Jo L ^ Jo Jo 



(37) 



If Po(^) ~ ^ as t> — * 0, Eq. (37) suggests a similar behavior for the steady state velocity 
distribution function, 

P 00 (v)^Av u i><1. (38) 

Since the concentration of clusters tends to the steady state limit = / °° Poo(v) dv, the 
exponent v must satisfy the inequality v > — 1. If one substitutes the assumed power law 
behaviors for Po(^) an d Poo(v) into Eq. (37), three possibility arise in the limit v — > 
which depend on the sign of v — ^ + 1. In the case where z/ > f — 1, the first exponential 
factor in Eq. (37) provides the dominant contribution. However, a simple calculation of 
the integral shows that v = f — 1- Similarly, for v < f — 1, one again finds v = f — 1- 
Only the last possibility, v = ^ — 1, appears to be self-consistent. Since v > —1, we obtain 
fj, > 0. Therefore starting from the assumption that the system reaches the steady state 
we have obtained that the exponent \i should be positive. This provides evidence for our 
conclusion that fj, = demarcates the scaling and steady state behaviors. 
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Notice also that for < /i «C 1, an asymptotic analysis of Eq. (37) gives the nu- 
merical prefactor in Eq. (38), A ~ This yields the estimate for the steady state 
concentration of clusters, 

(2a 

< \i < 1, (39) 




Since diverges as \x — > this indicates that at the critical value /z, ji = 0, the system is 
still evolving. 

Assuming the scaling form, obtained by the power counting analysis of Eq. (36) , let us 
examine the asymptotic behavior of the velocity distribution and the typical concentration. 
If we substitute the scaling assumptions 



P(v,t)~tP- a il>(V) with a = 



fi 



(40) 



into Eq. (36), we arrive at the following equation for the scaling function ip(V), 

,i r -, r v 

iP(V) = aV>* rfrexp — br 2 V^+ 2 - r(l - rf~ a / ^[V'(l - rf](V - V)dV 
Jo 2 J 

"(41) 

where t = t' ft and 6 = a/(/z + l)(/z + 2). 

Although we are unable to solve this nonlinear integral equation in general, we can 
obtain interesting information regarding the interesting borderline case of \x = 0. This case 
corresponds to a = 0, suggesting that the concentration grows slower than algebraically 
in time. However, the concentration of clusters is given by 



c(t) ~ t~ a J i()(V)dV, 



(42) 



and for \i = the integral diverges at the upper limit. To obtain the asymptotic behavior 
we consider Eq. (41) for the case n = 0, 



fi [ 1 rv 

ifj(V)=a rfrexp — br 2 V 2 - r(l - r) / ip[V'{l - t)](V - V)dV 
Jo 2 J 



(43) 



If we temporarily ignore the second term in the exponent, we find ip(V) ~ V~ x for V — > oo. 
If we then include the second term in the exponent and apply the previous asymptotic 
behavior of ip(V) ~ V^ -1 , we find tp{V) ~ (V logF) -1 . These estimates suggest the ansatz 



i/>(V) -CV- 1 (logV)~ x , 



(44) 
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for V >> 1. Upon substituting Eq. (44) into Eq. (43) one obtains the constants A = 1/2 
and C = l/y/2. With these values, the integral in the right-hand side of Eq. (42) diverges 
as ^2 log V, where V now denotes the maximum value of the scaled velocity. Since f3 = 1, 
this maximal velocity is proportional to t, which therefore suggests the logarithmic time 
dependence 

c(t) ~ y/2 logt. (45) 

In the complementary case of fi < 0, we have confirmed that the naive scaling ansatz is 
consistent with Eq. (42). 

VII. Conclusion 

We have introduced a simple ballistic aggregation process that mimics the kinetics of 
clustering in a single lane of traffic. Through direct probabilistic approaches, the analytical 
forms of the cluster velocity distribution and the joint mass- velocity distribution have been 
derived. For an initial velocity distribution of the form Pq(v) ~ as v — > 0, both the 
average velocity and the average cluster mass have power law time dependences with 
exponents that are rational functions of fi. This qualitative behavior is similar to that 
observed in the closely related ballistic annihilation process. We are also able to determine 
the asymptotic form of the joint mass-velocity distribution. 

Our model can also be analyzed in the cases of a spatially heterogeneous particle 
distribution and continuous input of particles. For the simple case of an initial one-sided 
spatial distribution, the system evolves towards a scaling distribution both in velocity and 
spatial variables. We have thus found that the total number of clusters in the initially 
empty half-line grows logarithmically with time for all initial velocity distributions. When 
there is a steady input of particle in an initially empty system, we have found that there 
is a transition between steady state behavior for fi > and transient behavior for fi < 
which is similar to that found when there is no in put. For the borderline case of the 
uniform distribution, fi = 0, we have found that the total concentration of clusters grows 
as \/logt. 

The irreversible traffic model introduced in this paper leads to ever-growing clusters. 
To describe traffic flows more realistically, several mechanisms to induce a steady state 
can be envisioned. For example, the input model can be generalized to incorporate a flux 
out of the system. Another realistic direction is to allow a faster car to pass a slower car 
at a rate which is some increasing function of the velocity difference of the two cars. This 
would allow a fast car to traverse a cluster car-by-car and ultimately regain its intrinsic 
velocity once the cluster is completely passed. It may prove interesting to examine the 
steady-state properties for this class of models. 
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Appendix A: Derivation of Eq. (24) 

In this appendix, we derive the asymptotic form of the cumulative velocity-mass dis- 
tribution which is valid for an arbitrary initial velocity distribution with a power-law small 
velocity tail. The starting point for calculation of Q m (v,t) is the formal expression of 
Eq. (23). By interchanging the order of the velocity and spatial integrations, the expres- 
sion can be rewritten as 

(.),„( c.l) Pie. I) ] [ / dxiexp(-Xi) dviP (vi). (Al) 

i=1 Jo Jv+{x^ VXi)/t 

Since the details of the initial distribution for large velocities do not change the form of 
the scaling solution, we may treat the more general power-law case by choosing a specific 
initial distribution of velocities whose form is convenient for performing the integration 
over the velocity in the right-hand side of Eq. (Al). Since the velocity distribution near 
v = has a power-law tail, the most convenient initial distribution is 

P (v) = av 11 exp {-av^ 1 / {n + 1)) . (A2) 

For this initial distribution, the integration over the velocity variables is immediate and 
one finds, 

m — 1 «oo r jU+1 

Q m (v,t) = P(v,t) Y[ / dxiexp(-Xi) exp -a(v+(x 1 -\ h /(a* + 1) 

i=i J° L 

(A3) 

To evaluate the multiple spatial integral, we define f(z) = exp (-az» +1 /(iJ,+ 1)) 
and expand f(z) to first order about the point z = v and then exploit a number of 
simplifications associated with performing the integrals over the factors e~ Xi . This gives 

m-l „oo 

Qm(v,t)=T] / dx z e~ Xi f (v + e(xi H + x z )) 

i=i Jo 

m-l „oo 

~ TT / rfx l e-^/(z;) + (e(x 1 + --- + x l ))/ / ^), (A4) 
i=i Jo 

m — 1 

w H /(^ + «)+0(e 2 ), 



where e = 1/t. By substituting the explicit functional form, f(z) = exp(—az^ +1 /(fi + 1)), 
we have 



m — 1 



Q m (v, t) = P(v, t) 11 exp (-a(v + i/ty +1 /(pi + 1)) . (A5) 



i=l 
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Finally, the product of the exponential factors is written as an exponent of a sum. In 
the limit of large m, this sum is equivalent to the integral a f™ 1 dy(v + y/t)^ +1 /(/z + 1). 
Evaluating this integral, the asymptotic form of the cumulant mass-velocity density is 
obtained as 

Q m (v, t) ~ av» exp [-b{vt p + mr a y +2 ] . (A6) 

In evaluating this asymptotic expression, the exponential factor of P(v, t) cancels the factor 
that emerges from lower limit of the integration over y. In Eq. (24), the above expression 
is written as a scaling function. 

Appendix B: Analytical Solution for the Exponential 
Initial Velocity Distribution 

We outline here the explicit analytical solution for Q m (v,t) for the exponential initial 
velocity distribution, Pq(v) = e~ v , which corresponds to the special case \x = and a = 1. 
For this exponential distribution, the integration of Po(vi) over the velocity variables Vi, 
according to Eq. (Al), equals exp[— (v + {x\ + • • • + Xi)/t)\. Hence, we obtain 



m— 1 „oo r / _ • \ 

Q m (M) = P(M)e- (m - 1)u J] ^ dx i exp\-x i (l + — . (Bl) 



Upon integration over the space variables, the following exact expression is found for 
Q m (v,t), 

Q m (v, t) = P(v, t) e -("-i>«t"-i ^ ± 1} , (B2) 

1 (t + m) 

where T(z) = / °° x z ~ x exp(— x) is the Euler gamma function and the velocity distribution 
obtained from Eq. (14) is P(v,t) = exp[— v — t(e~ v — 1 + v)]. 

The exact forms for the joint mass-velocity distribution and for the mass distribution 
can be evaluated from Eq. (B2) by taking the appropriate limits m — > oo and v — > 0. The 
resulting asymptotic expressions are identical with the expressions of Eqs. (25) and (26) 
respectively. 
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Figure Captions 

Figure 1. Schematic illustration of the irreversible traffic model. A faster car overtakes 
a slower car and after the encounter, the faster car assumes the velocity of the 
slower car. 

Figure 2. Illustration of the initial configuration of a possible three car cluster. 
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